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Abstract 


A detailed theoretical and numerical investigation of the infinitesimal single-crystal gradient- 
plasticity and grain-boundary theory of Gurtin (2008) “A theory of grain boundaries that accounts 
automatically for grain misorientation and grain-boundary orientation”. Journal of the Mechan¬ 
ics and Physics of Solids 56 (2), 640-662, is performed. The governing equations and flow laws 
are recast in variational form. The associated incremental problem is formulated in minimiza¬ 
tion form and provides the basis for the subsequent finite element formulation. Various choices 
of the kinematic measure used to characterize the ability of the grain boundary to impede the 
flow of dislocations are compared. An alternative measure is also suggested. A series of three- 
dimensional numerical examples serve to elucidate the theory. 
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1. Introduction 


The miniaturisation of mechanical components composed of crystalline material requires a 
continuum theory that accounts for the role of the grain boundary and for size-dependent effects. 
The grain-boundary model should incorporate both the misorientation in the crystal lattice be¬ 
tween adjacent grains, and the orientation of the grain boundary relative to the crystal lattice 
of the adjacent grains. Classical theories of plasticity are unable to describe the well-known 
size-dependent response exhibited by crystalline material at the micro- and nanometre scale. 
Numerous extended (gradient and non-local) continuum theories of single-crystal plasticity have 
been presented in the last two decades to circumvent these limitations. The thermodynamically 
consistent gradient theory of Gurtin and co-workers and related works (see e.g. Gurtin, 2002, 
2006, 2008a; Gurtin and Needleman, 2005) have received particular attention. A variational for¬ 
mulation of the Gurtin (2002) framework has been developed in Reddy (2011a,b). In Gurtin 
(2002) the defect part of the free energy is parametrized in terms of the (bulk) Burgers tensor, 
a rigorously defined and physically meaningful measure of the (local) Burgers vector and hence 
the lattice mismatch (see e.g. Nye, 1953). The form of the defect energy was modified by Gurtin 
(2008a) to account for a continuous distribution of geometrically necessary dislocations (GNDs). 
The recent work of Gurtin and Reddy (2014) uses a scalar measure of the accumulated slips as 
the basis for the hardening relation, which takes account of both self- and latent-hardening. Fur¬ 
thermore, the resulting initial boundary-value problem is placed in a variational setting in the 
form of a global variational inequality. Ertiirk et al. (2009) show how the theory of Gurtin et al. 
can be related to the more physically motivated theories due to Bayley et al. (2006); Evers et al. 
(2004a,b, 2002). 

The gradient theory of Gurtin (2008a) provides a basis to account for the role of the grain 
boundary (see Gurtin, 2008b). Neumann and Dirichlet-type boundary conditions on the slip 
and the flux of the vectorial microforce (i.e. the microscopic traction), respectively, can be pre¬ 
scribed and are often assumed homogeneous. The homogeneous Dirichlet condition, known as 
the micro-hard boundary condition, has been widely used to account for the grain boundary or 
an interface (see e.g. Ekh et al., 2007; Evers et al., 2004b; Kuroda and Tvergaard, 2008; Lele and 
Anand, 2008; Ohno and Okumura, 2007). Clearly this boundary condition ignores the complex 
geometric structures in the vicinity of the grain boundary. 

Central to the theory of Gurtin (2008b) is the introduction of the grain-boundary Burgers 
tensor to parametrize the grain-boundary free energy. The grain-boundary Burgers tensor is 
obtained from the mismatch in the plastic part of the displacement gradient around a circuit 
centred on the grain boundary and contains information on both the misorientation in the crystal 
lattice between adjacent grains, and the orientation of the grain boundary. Eurthermore, the 
definition of grain-boundary Burgers tensor is consistent with that of the Burgers tensor in the 
bulk. 

The grain-boundary Burgers tensor can be expressed in terms of the intra- and inter-grain 
interaction moduli. The inter-grain interaction moduli account for mismatch in the slip sys¬ 
tems adjacent to the grain boundary and the orientation of the grain boundary. They provide a 
physically meaningful characterisation of the interaction of neighbouring slip systems with the 
extremes described as non-interactive and maximally interactive. 

Recently, van Beers et al. (2013) have proposed and numerically implemented a theory simi¬ 
lar to that of Gurtin (2008b) for incorporating grain boundaries into the gradient crystal plasticity 
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theory of Evers et al. (2004b). A key feature of the theory is the use of a geometrically-motivated 
vectorial measure to parametrize the grain-boundary free energy. We will show in this work that 
under planar conditions (as investigated numerically by van Beers et al. (2013)) the models of 
Gurtin (2008b) and van Beers et al. (2013) produce identical interaction moduli. In similar work, 
Ozdemir and Yalginkaya (2014) implemented the grain-boundary theory of Gurtin (2008b). A 
series of finite element simulations of planar bi-crystals (single and double slip) illustrated fea¬ 
tures of the grain-boundary model. 

Gudmundson (2004) and Fredriksson and Gudmundson (2005a,b) propose an interface the¬ 
ory in which both the interface moment traction and the plastic slip can be discontinuous at 
the grain boundary. The model introduces an interfacial free energy that depends on the plastic 
strain on both sides of the interface. Critically however, the model does not account directly 
for the mismatch in the adjacent grains or the orientation of the grain boundary. Related works 
include those by Aifantis and Willis (2005, 2006). Ekh et al. (2011) propose a “micro-flexible” 
grain boundary which provides a degree of resistance to plastic flow dependent upon the mis- 
orientation of the adjacent grains (see Bargmann and Ekh, 2013, for an extension that accounts 
for thermal effects). They do not consider the orientation of the grain boundary. Wulfinghoff 
et al. (2013) account for grain boundaries within a gradient-plasticity theory by postulating a 
grain-boundary yield condition and flow rule. The theory does not account for the mismatch in 
the adjacent grains or the orientation of the grain boundary. Voyiadjis et al. (2014) developed 
and numerically implemented a gradient-plasticity model for the polycrystalline problem which 
accounts for the role of grain boundaries via the mismatch in the accumulated plastic strain. 

Recasting the problem of single-crystal gradient plasticity as a variational formulation makes 
it amenable to analysis (see Reddy, 201 lb). The variational formulation does not have an associ¬ 
ated minimization problem, but the corresponding time-discrete incremental problem does. The 
variational formulation developed in Reddy (201 lb) is extended here to include the grain bound¬ 
ary and a viscoplastic flow law. The associated incremental minimization problem is shown to be 
equivalent to the time-discrete variational problem and provides the point of departure for the nu¬ 
merical implementation within the finite element method. The software AceGen (Korelc, 2002) 
is used to describe the finite element interpolation, and to compute the residual and (algorithmi¬ 
cally consistent) tangent contributions directly from the prescribed functional associated with the 
incremental minimization problem, using automatic differentiation at the level of the quadrature 
point. This approach ensures quadratic convergence of the algorithm and greatly simplifies the 
implementation. Details of the numerical implementation are given. 

Gurtin (2008b) proposes two thermodynamically admissible plastic flow relations for the 
grain boundary (denoted Gurtin I and II). The flow relations define the structure of the dissipative 
microscopic stress in the grain boundary microscopic force balance. The flux of dislocations 
from the grains drives the microscopic force balance. In the first proposal, the grain boundary 
Burgers tensor is used to parametrize the flow relations, while in the second it is the slip. The 
first approach accounts for the interaction of slip systems adjacent to the grain boundary. This 
approach also allows for a recombination of the plastic distortion contributions from adjacent 
sides via the definition of the grain boundary Burgers tensor. The second approach does not 
directly account for the structure of the adjacent grains or the orientation of the grain boundary 
in the plastic flow relation. Both approaches account for the geometric structure of the adjacent 
grains and the grain boundary via the flux terms from the grains. 
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A series of three-dimensional numerical examples elucidate the grain-boundary theory. The 
examples demonstrate single slip in a bi-crystal and multi-slip in a polycrystal where each of 
the 27 grains is a face-centered-cubic crystal structure. The polycrystal example in particular 
demonstrates various features of the Gurtin (2008b) theory that are not obvious from the theory 
or the single slip examples. The Gurtin I model for the plastic flow relation is unable to capture 
the widely-used micro-hard condition in multi-slip problems, even when using an artificially high 
value for grain boundary slip resistance. The Gurtin II model can capture the range of responses 
between the micro-free and micro-hard conditions. Motivated by a plastic flow relation that 
accounts for the structure of the grain boundary and captures the micro-hard and micro-free limits 
and the range between, a modified measure of the grain-boundary Burgers tensor is analysed and 
implemented. In the modified formulation, the micro-hard limit is recovered for large-angle grain 
boundaries, and the micro-free for perfectly aligned crystal structures on either side of the grain 
boundary. 

The structure of this work is as follows. The kinematics of the gradient crystal plasticity 
formulation in the bulk and on the grain boundary are reviewed in Sec. 2. The kinematic measures 
of the mismatch at the grain boundary proposed in Gurtin (2008b) and van Beers et al. (2013) 
are then compared. Particular attention is paid to the inter-grain interaction moduli. The kinetics 
of the problem and the various restrictions to the theory assumed are presented in Sec. 3. The 
governing relations and the plastic flow law are presented in Sec. 4 and 5. An alternative measure 
of the kinematic mismatch at the grain boundary is given. The variational formulation of the 
problem and the associated incremental formulation are developed in Sec. 6. This is followed by 
details of the numerical implementation within the finite element framework. The finite element 
model is then used to simulate a series of representative numerical examples in Sec. 8. Finally, 
conclusions are made and various extensions proposed. 

Notation and basic relations 

Direct notation is adopted throughout. Occasional use is made of index notation, the sum¬ 
mation convention for repeated indices being implied. When the repeated indices are lower¬ 
case italic letters, the summation is over the range (1,2,3). Upper-case italic indices can re¬ 
fer to arbitrary adjacent grains (A, B}. The summation convention is not employed for grains. 
The scalar product of two vectors a and b is denoted a ■ b - [a],[ft],. The scalar product 
of two second-order tensors A and B is denoted A : B - [A]y[B]y. The composition of 
two second-order tensors A and B, denoted AB, is a second-order tensor with components 
[AB]ij = [A]im[B]^j. The tensor product of two vectors a and ft is a second-order tensor 
D - a ® b with [Djy = [«],[ft];- The action of a second-order tensor A on a vector ft is a 
vector with components [a], = [A],„[ft]m. The curl of a second-order tensor A is a second- 
order tensor with components [curlA],^ = eirsdAjJdXr, where e is the third-order permutation 
tensor. An arbitrary quantity in the bulk is denoted {•) and analogously {•) denotes an arbitrary 
quantity on the grain boundary. Any array associated with the set of N slip systems is denoted 
7 ■- ^ 7^ ■ ■ ■ >7^1 • Summation over the slip systems will be abbreviated by Yja- Index nota¬ 

tion is not employed for summations over slip systems. The unit basis vectors in the Cartesian 
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(standard-orthonormal) basis are {ei, € 2 , e^]. The following identities are employed widely: 


[wX]/y .— ^ikj^k ■> 

-[nx][nx] - I - n® n, 

[nx]b - nxb, 

where I is the second-order identity tensor. 

2. Kinematics 

The kinematic basis of the gradient-plasticity and grain-boundary theory of Gurtin (2008b) 
is first summarised. A comparison of aspects of the theory with that of van Beers et al. (2013) is 
then made. The theory presented is based on the assumption of infinitesimal deformations. The 
initial configuration is thus assumed to be geometrically representative for all time. 

2.1. Bulk 

Consider a continuum body whose placement 'V at time t = 0 is shown in Fig. 1. A typical 
material point is identified by the position vector x. The displacement of the material point is 
denoted by u(x, t). The displacement gradient // := Vm is decomposed (locally) into elastic and 
plastic parts and according to 


The elastic displacement gradient accounts for recoverable elastic lattice stretching, while IP 
quantifies the plastic distortion due to slip on the predefined slip planes. The elastic strain is 
given by the symmetric part of the elastic displacement gradient 

\ -h . 

The flow of dislocations through the crystal lattice is described kinematically via the assump¬ 
tion that the plastic distortion tensor can be expressed in terms of the slip 7 “ on the individual 
prescribed slip systems a — 1 , 2 ,..., A as 

TfP = ^ ® ^ r“S“. ( 1 ) 

a a 

The slip direction and slip plane normal of slip system a are denoted and /m“, respectively, 
where s“ ■ nf =0 and |x“| = \nf\ - 1. The vector /“ is defined by /“ nf x Hence 
{nf,P] constitute a local orthonormal basis. The Schmid tensor is defined by S“ = s® ® m". 
The Burgers tensor 

G = curlTfP = ^ [Vr“ X nf] ® , (2) 
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Figure 1: A continuum body ‘V, composed of grains separated by a two-sided grain boundary Q. 

quantifies the crystal distortion due to dislocations. It is obtained from the boundary integral of 
IP over an infinitesimal closed circuit in the bulk. The vector G^e gives the Burgers vector, per 
unit area, for a closed circuit on a plane with unit normal e. 

Following Gurtin (2008a), the constitutive theory at the microscopic scale accounts for a 
continuous distribution of GNDs. The dislocations are either of edge or screw type and are 
characterised in terms of their Burgers and line directions as follows: 

• edge dislocation: Burgers direction s" and line direction F-, 

• screw dislocations: Burgers direction s“ and line direction s". 

The density of the edge and screw dislocations per unit area are related to the slip gradient as 



Po = 


and 


(3) 


where b is the length of the Burgers vector. 

2.2. Grain boundary 

The body 'V is composed of grains separated by grain boundaries as depicted in Fig. 1. The 
boundary of a grain interior to the domain is denoted by The set of all grains is denoted 

by V = {'Va}, with 'Va an arbitrary grain. The grain boundary is represented as a two-sided 
interface, with common sides 0a and Qg. The set of all grain boundaries is denoted by G = {0c)^ 
where C runs from one to the number of grain boundaries. The grains on either side of the grain 
boundary are correspondingly denoted by 'Va and 'Vb- The outward unit normals to 0a and 0b 
are denoted Ha and n^, respectively. The normal to the grain boundary is defined by n := Ha- 
The displacement field across the grain boundary is assumed to be continuous (see Gurtin and 
Anand, 2008, for a model of grain boundaries that slip and separate); that is. 


Ih]1 := u\g, - U\g^ = 0 


and 


u u\g. 


(4) 


6 





The plastic slip and its gradient are however generally discontinuous across the grain boundary. 
Central to the grain-boundary theory of Gurtin is the grain-boundary Burgers tensor 

obtained by integrating the plastic distortion //p along an infinitesimal circuit centred on the grain 
boundary (cf. Eq. (2)). The grain-boundary Burgers tensor can be expressed in terms of the slip 
on either side of the grain boundary, the orientation of the slip systems and the grain-boundary 
orientation as 


a a 

where 


Wj = S^[nx] = ® [/m“ X n] 


for 7 = A, B, 


is the Schmid orientation tensor for system a in grain J. 

The Burgers vector production within the grain boundary is given by 


G = H^nx] = [fsK - G“. 


(5) 


The contraction of the Schmid orientation tensor with the grain-boundary Burgers tensor 


motivates the definition of the slip-interaction moduli := : N^, that is, 

G = y] [r“ [n^ Nl]] 


where 



Furthermore, 

C“j - [«“ • [nfj X n] ■ [ttiPj X n]. (7) 

The slip interaction moduli are classified as: 

• Cf,: intra-grain interaction moduli for grain A; 

• Cf: inter-grain interaction moduli between grains A and B; 

• cf: intra-grain interaction moduli for grain B. 

The inter-grain interaction moduli for an arbitrary pair of slip systems and grain-boundary ori¬ 
entation thus have the range -1 < Cf < 1. The inter-grain interaction moduli characterise the 
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interaction between slip systems in adjacent grains separated by a grain boundary. The systems 
are non-interactive if = 0 and maximally interactive if = 1. These definitions imply 
that a pair of slip systems in adjacent grains are (see Eq. (7)): 

• non-interactive if and only if either: 

■ 5^ = 0 or [»i^ X n] ■ [»4 x n] = 0; 


• maximally interactive if and only if all of the following conditions are satisfied: 

si - ±s1 , m1-n-0, mPg-n-Q. 


A key measure in the Gurtin (2008b) model used to parametrize the grain-boundary free 
energy is the norm of the grain-boundary Burgers tensor: 


IGP 


G:G = £g“| 


( 8 ) 


which is a function of the adjacent slip system structures and the grain-boundary orientation. 


2.3. Comparison of aspects of the Gurtin and van Beers et al. models 

The kinematic measures of the mismatch at the grain boundary employed in the models of 
van Beers et al. (2013) and Gurtin (2008b) are now compared. The model of van Beers et al. 
(2013) is stated using notation similar to Gurtin (2008b) to aid comparison. 

Van Beers et al. (2013) define the geometrically-motivated grain-boundary normal slip com¬ 
ponents (see Ertiirk et al., 2009; Kuroda and Tvergaard, 2008) by 

Ta K - ^a] ■ « and q% := 7 “ [-s“ + I1]-n. 

The net defect vector (modelling parameters are set to unity) on the grain boundary is defined by 

g - > 


which can be restated in a form conducive for comparison as 

^ ® - Sb ® Sfi] « -Ta [*a ® ^a - *a ® *a] « 

^ III" - X 

where |g| is used to parametrize the grain-boundary free energy (cf. Eq. (8)). The intra- and 
inter-grain interaction moduli in the van Beers et al. (2013) model are denoted by . 






+rArl' 


■^A 



A purported difference between the models of Gurtin (2008b) and van Beers et al. (2013) is in 
the characterisation of the inter-grain interaction moduli, and The relationship between 
G and ~g can be seen directly by assuming a planar problem in which /“ = —e^ and the slip plane 
and grain-boundary normal are in the e\-e 2 plane. Screw dislocations cannot be accommodated 
in such a model as /“ ■ n = 0. Contracting the grain-boundary Burgers tensor (using Eq. (9.8) in 
Gurtin, 2008b) with -e^ gives the net defect vector ~g\ 


-Ge^i — 


X l-ysK ■ ® ■ msi ® ® 11] + yiui ■ msi ® ^3 


-Yj^ylK-nK-ylK-nK] 


g- 


Hence the general theory of Gurtin (2008b) reduces to that of van Beers et al. (2013) under planar 
restrictions. Furthermore, 

^°AB = ■”][*!■ w] (planar problems). 

Consider the example of single slip in the bi-crystal shown in Fig. 2. The slip system in grain 
A is fixed. The slip system grain B is initially the same as A and is then rotated by an angle as 
around the ea-axis. The grain-boundary normal is initially at the value shown in Fig. 2 and then 
rotated by an angle ag around the e 3 -axis. Consider the planar case where the slip directions, 
slip plane normal, and the grain-boundary normal are constrained to lie in the ei-e 2 plane. As 
expected, the two theories produce identical inter-grain interaction moduli. 

The geometric arguments underpinning the vectorial quantity g are insufficient when the 
problem becomes three dimensional, and the two theories differ. The grain-boundary Burgers 
tensor G, however, is a fundamental measure that is well defined in three dimensions. It is 
clear from Fig. 2 that, as suggested by Gurtin (2008b), the inter-grain interaction moduli vary 
continuously between two extremes -1 < Cab ^ 1 and correspond to physically intuitive notions 
of non-interactive and maximally interactive slip systems. 

The computed intra-grain interaction moduli for the planar problem are identical. The intra¬ 
grain interaction moduli differ for non-planar problems. 


3. Kinetics 

The symmetric stress tensor in the bulk is denoted by T. Scalar and vector microscopic 
forces, denoted by n‘‘ and respectively, are postulated as conjugates to the slip rates and their 
spatial gradients (Gurtin, 2000, 2002). The vectorial microscopic stresses exist only in the bulk. 
The resolved shear stress on slip system a is denoted by t" T : S“. The power conjugate 
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Figure 2: Comparison of the inter-grain interaction moduli in the models of Gurtin (2008b) and van Beers et al. (2013) 
for a range of ag and ag [rad]. For the choice w = [0 1 0] the problem is planar. 


pairings are as follows; 




in^ 

(macroscopic stress). 

if 

in^ 

(scalar microscopic force). 

t> 

X 

in^ 

(vector microscopic force), 

—O' • at —o • tk 

^ ^ Jb 

on Qa,Qb 

(scalar microscopic forces). 


The scalar microforces on either side of the grain boundary are denoted by ;r“. 

3.1. Bulk 

The free energy in the bulk T* is decoupled into an elastic (macroscopic) and a defect (micro¬ 
scopic) part, denoted by 'T® and 'T'^ respectively, given by* 

T* = -hT'‘*(p) . 


*For further examples of defect energies and the inclusion of classical hardening in the free energy see Reddy (201 lb). 
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The defect energy is chosen to be the following function of the dislocation densities p: 


a 


where the material parameters are given by (see e.g. Ertiirk et ah, 2009; Evers et al., 2004b) 


^ 8 [l-v] 


and 




4 


and R is the back stress cut-off radius and v is the Poisson’s ratio. The dependence of on the 
gradient of the slip is via relation (3). 

The elastic free energy is chosen to be quadratic, for the sake of simplicity. The fourth-order 
isotropic elasticity tensor in the bulk is given by 

— 46ij6kl + P ? 


where A and p are the Lame constants. 

The standard Coleman-Noll procedure gives the stress as the conjugate kinetic quantity to 
the elastic strain: 





(9) 


The scalar energetic defect forces associated with edge and screw dislocations are respec¬ 
tively defined by 




d'¥‘'(p) 

'W~ 


and 


fa _ 

/g ” 


5'pd(p) 

dp^o 


The energetic vectorial microstress is defined by 






( 10 ) 


and the reduced dissipation inequality in the bulk follows as 




.O' I „o:‘a 

7T y 


> 0 , 


where is the dissipative vectorial microstress. This general framework is restricted (to elu¬ 
cidate the role of the grain boundary) to purely energetic vectorial microstresses and purely 
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dissipative scalar microforces: 


' ^ ” ^en 

_. „a _ „a 

^ ^ = ^dis 

^ ;r“r“ > 0 . (11) 

a 

3.2. Grain boundary 

Following the general theory of Gurtin (2008b), we consider both energetic and dissipative 
structures in the grain boundary. The influence of these two structures was numerically investi¬ 
gated by van Beers et al. (2013). 

The free energy of the grain boundary is parametrized in terms of the grain-boundary Burgers 
tensor, that is T* = T'(G). The grain-boundary energetic stress is defined by 


n‘L = 0 


M 


5'T(G) 

dG 


thus 


T' = M : G. 

Following Gurtin (2008b), the grain-boundary free energy is assumed to be given by the quadratic 
function 


T'(G) = i/l|G|2, 


which yields 


M = TG, 

where d > 0 is a constant constitutive parameter. The rate of change of the grain-boundary free 
energy can be expressed in terms of the slip rates on either side of the grain boundary, using 
Eq. (5), as follows: 


T' = M:G = dG:G 


= ;^[[dNg:G]r2-[^N^GK] 

a 



a 


where the energetic internal microforces are defined by : = d : G and 7r“ : = -d : G. 

Following Gurtin (2008b); Gurtin and Needleman (2005), the reduced dissipation on the 
grain boundary is expressed in terms of the rate of change of the grain-boundary Burgers tensor 
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G and its conjugate kinetic quantity K: 


K:G>0, (12) 


where the dissipative internal microscopic forces on the grain boundary 7rA,dis and ttb, dis are de¬ 
fined by the relations 

and (13) 

4. Governing equations 

The equations governing the response of the bulk and the grain boundary are now sum¬ 
marised. 

4.1. Bulk 

A balance of macroscopic and microscopic forces in the bulk yields 

divT = 0 in^ and t*{n)-Tn on&Vti, (14) 

div^“ H-t“ - ;r“ = 0 in^^ and E“*(n) = ■ n on 5^^;-. (15) 

Eq. (14) is the standard equilibrium equation in the absence of body and inertial forces, and t* is 
the prescribed Cauchy traction on the Neumann part of the boundary Dirichlet boundary 

conditions on the displacement u are prescribed on where U and d^N H 

&V 0 - 0. Furthermore, the boundary &V is subdivided into complementary parts and 
d^H- The standard boundary condition on the micro-free part of the boundary d'Vp is that 
the scalar microscopic traction = 0 while on the micro-hard part of the boundary d^H 
homogeneous conditions on the slip are prescribed. For additional details on the microscopic 
boundary conditions see Gurtin and Needleman (2005). The macroscopic and microscopic force 
balances are coupled via the dependence of the resolved shear stress on the macroscopic stress 
T. 

4.2. Grain boundary 

The balance of macroscopic and microscopic forces across the grain boundary yields 

irp = 0 onG (16) 

= 0 on^A and ■ n = 0 on Gb ■ (17) 

Eq. (16) is the standard traction continuity condition for an interface. The microforce balance 
(17) on the grain boundary states that the internal microforce on either side of the grain boundary 
acts in response to the flux of the vectorial microforce from the grain. The gradient-plasticity 
formulation adopted in the bulk allows meaningful balance equations for the grain boundary to 
be constructed in a consistent manner. A micro-hard or micro-free boundary condition could 
be applied on the grain boundary. However all information about the geometry of neighbouring 
crystal structures and grain-boundary orientation would be lost. 
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A,dis ^ o,dis 



5. Plastic flow relations 


In order to complete the theory, the reduced dissipation inequalities in the various parts of the 
body (11, 12) need to be satisfied in a thermodynamically consistent manner. The plastic flow 
relations are then obtained by postulating forms for the dissipative internal variables. 


5.1. Bulk 

The yield function /(7r“) defines the region of admissible microscopic forces on the a* slip 
surface. The yield function and the flow law for the plastic slip are defined by 

/(7r“) = |;t“|-5 <0, (18) 

dfin") 

r = =^“sgn;r“, (19) 

a;r“ 

where 5 > 0 is the constant slip resistance and /i“ > 0 is a scalar multiplier, together with the 
Kuhn-Tucker complementarity conditions 

/(O<0, /l“>0, d“/(7r") = 0. (20) 


Under conditions of plastic flow /(7r“) s 0 and the flow rule (19) can be inverted to obtain 


n" = Ssgny" = 


dD' 


eft 


(9y“ 


( 21 ) 


where the effective dissipation function is given by 51'/'l. More generally, the elastic region, 
the flow rule and the complementarity conditions can be expressed in the alternative, equivalent 
form 


z):(f(r“)>Dyt“) + 7r«[r“-r“]- 

This formulation, which has received a detailed treatment in Han and Reddy (2013), is less 
popular as a basis for computational treatments of the classical problem than (18)-(20), which 
uses the flow law in its traditional form of the normality law with the Kuhn-Tucker conditions. 
Nevertheless this formulation is particularly well suited to problems such as gradient plasticity 
in which the plastic slip is a primary unknown. 

The rate-independent theory presents various numerical challenges due to the indeterminacy 
of plastic slip. To circumvent these problems, a regularized effective dissipation function is 
proposed of the form 


0“ = - 

m+\ 


to 


im+l 


ro, 


where fo is the reference value for the slip rate and m > 0 is the rate sensitivity. The scalar 
microforce follows as (cf. Eq. (21)): 


TT 


a 


dD". 

VIS 

dy" 





sgnt“ • 


( 22 ) 
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5.2. Grain boundary 

The grain boundary impedes the flow of dislocations from the adjoining grains. The geomet¬ 
rical structure of the grain boundary quantifies the amount of impedance. Should a measure of 
the microforce on the grain boundary reach a sufficient level, the grain boundary can yield and 
allow for transmission of dislocations. 

Gurtin (2008b) proposes two forms for the flow relation on the grain boundary, both of which 
satisfy the reduced dissipation inequality on the grain boundary (12). The first is a flow relation 
explicitly parametrized by G (here referred to as the Gurtin I model). The flow relation therefore 
accounts for the interaction of the slip systems in adjacent grains and the grain boundary orienta¬ 
tion. The second relation is parametrized in terms of the slip on either side of the grain boundary 
(here referred to as the Gurtin II model). 

The Gurtin I model allows for the interaction of slips from either side of the grain boundary 
via the definition (5) of G. It is thus possible to have a G s 0, and hence zero dissipation, for 
non-zero values of the slip rate y/. The Gurtin II model excludes the possibility of non-zero slip 

resulting in zero dissipation. A modified definition for G is given in Sec. 5.3 that, when used to 
parametrize the flow rule, ensures that G s 0 if and only if ji = 0. The response of the three 
models are explored via a series of numerical examples in Sec. 8. 


5.2.1. Plastic flow relation in terms ofG (Gurtin I model) 

A regularized dissipation function for the grain boundary is chosen as 




m + I 


to 


m+1 


Go, 


(23) 


where 5 > 0 is the slip resistance. Go is the reference value for the rate of the grain-boundary 
Burgers tensor and m > 0 is the rate sensitivity parameter. The flow relation for K follows as 


dG 


1 ^ 

^0 


G 

|G| 


(24) 


It follows from Eq. (13) that 




to 


: G 


and 




1^ 

^0 


: G 

|G| 


The above relations can be expressed in terms of the slip rates using Eq. (5). 


5.2.2. Plastic flow relation in terms off (Gurtin II model) 

A regularized dissipation function for side I of the grain boundary is chosen as 



S 

m H- 1 



To 


To- 
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where jq is the reference value for the rate of the grain-boundary slip, and the remaining pa¬ 
rameters are defined as in Eq. (23). The regularized dissipation function for the grain boundary 
follows as 


^vis ~ [^vis,A • 

a 

The dissipative internal microscopic forces on either side of the grain boundary (see Eq. (17)) 
are thus obtained as 


\r,\ ■ 

5.3. A modified form for the flow relation on the grain boundary 

Consider the Gurtin I model. It is possible to construct situations such that 

^ = Z = 0 (25) 

a a 

for non-zero values of on both sides of the grain boundary (see Appendix A for an illustrative 
example). Such situations are illustrated in Sec. 8.3 for a three-dimensional example with a face- 
centered-cubic crystal structure. This situation renders the expression (12) for the corresponding 
dissipation meaningless. It should be emphasised that this is not a deficiency of the Gurtin I 
model. It could be argued that the recombination of slip in the grain boundary accounts for 
physical phenomena such as annihilation. As demonstrated in Sec. 8, the model is not, however, 
capable of reproducing the behaviour of both the micro-hard and micro-free boundary conditions. 
With this as motivation we revisit the definition of the grain-boundary Burgers tensor with a view 
to constructing an alternative measure of “defect” density across the boundary. 

The first observation is that the Burgers tensor in the bulk G is obtained from the plastic 
deformation tensor IP, which by Eq. (1) involves a sum over all slip systems: the derivation of 
G involves taking a circuit of //p and then using Stokes’ Theorem. The derivation of the grain¬ 
boundary equivalent involves taking an infinitesimal circuit that includes a section of the grain 
boundary. In contrast to the situation in the bulk, the slip systems on either side of the boundary 
bear no relation to each other. In particular, there is no relationship between the numbering 
adopted on either side. So instead of considering the jump in IP, that is, 

a 

we consider the differences between the individual components that make up IP. 

By analogy with the definition of G“ consider the difference 



Here is a flag whose components take the values 1 or 0, in this way allowing a selection of 
particular combinations of slip systems for the two adjacent grains. The square root is introduced 


‘/,dis 






It? 


Tn 
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in order to arrive at a convenient formulation for the magnitude later. Then 






( 26 ) 


We note also from Eq. (25) that 




aa 
6 ’ 


in which the subscript 6 denotes the choice = (5 q,^. 

Now define a modified magnitude Er of G by 

This may be compared with the conventional magnitude which is given by 

|2 


|G|^- 


-laa 

’d 


Z 6 : 

a 


Thus we see that the conventional magnitude is the absolute value of a sum while the modification 
is the sum of weighted absolute values. 

We have 






2 

^<5 • 




Note also that for the special case in which = Sap, with Sap being the Kronecker delta, 


t^2 \ ' l/^orO'E 

= 2^ IG^ I • 

a 


Thus 


|G|2<r2 


5.4. The dissipation function 

The modified definition of the magnitude of G alters the form of the dissipation inequality and 
flow rule on the grain boundary. In what follows we drop the subscript R on for convenience. 
Eurthermore we consider only dissipative microforces on the grain boundary. 

The reduced dissipation inequality on the grain boundary (12) is restated as 

^ : G > 0. 

a p 
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Following a similar procedure that lead to the definition of K, we write 

+ : G“^ > 0, 

a p a,fi 

where has to be defined. Expand this expression to obtain 

^ K‘^13 ; 

aJ3 

= Z[Z 



where ^ and^ are the modified scalar microforces on adjacent sides of the grain boundary. 
For the rate-independent case the modified yield function / is defined by 

afi 

Then the normality relation is 

QRo^j} ■ 

The corresponding dissipation potential is 

£)(G“^) = 5f = 5 ^|G I 

afi 

SO that when flow occurs 

_ dD{G‘^P) 

dG"P f 

Remark. The modified and original formulations are identical for the case of single slip with 
Ra/s _ 2 -pjjg choice made here for in multi-slip problems is as follows. For each slip system 
a in grain A, the most interactive slip system in grain B is identified. The most interactive slip 
system in grain B is the slip system most likely to receive dislocations transmitted from slip 
system a in grain A, should plastic flow occur. The process is repeated for grain B. The simple 
example in the following section provides further details. Many other choices are possible; the 
choice made here is based on the concept of inter-grain interaction defined in the original theory 
of Gurtin (2008b) and ensures that the dissipation at the grain boundary is zero if and only if the 
slip on all adjacent slip systems is zero. The modified model contains features of both the Gurtin 
(2008b) I and II models. 
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5.5. An illustrative example 

The procedure to select the most interactive slip system in adjacent grains, A'a and A'b, 
separated by a grain boundary Q and hence compute 7?“^ is illustrated using the simple example 
problem shown in Fig. 3. The matrix of inter-grain interaction moduli Cab is computed using 
Eq. (7) as 


Cab = 


0.5 

1 


0 

0.5 


and 



1 

0.5 


Recall that from the symmetry relation (6) Cab - C^^. The matrix Cab, with components 
[Cab]“^, contains the interaction coefficients for all slip systems a in A'a with all other slip 
systems /? in ^b- 

Consider a slip system a in A'a- The slip system in A'b that interacts the most with a in A'a is 
given by the choice of /? that gives the maximum value of [Cab]“^- For example, slip system 1 in 
A'b interacts the most with slip system 2 in A'a- Consider now slip system a in 'Vb- Slip system 
2 in A^b interacts the most with slip system 2 in 'Va- The corresponding interaction matrices Rjj 
have a value of 1 in the ojS slot if slip system J3 in'Vj interacts the most with slip system a in A':, 
and 0 otherwise. The interaction matrices for each grain are 


^AB - 


0 

0 


and 


^BA - 


0 

0 


1 

1 


The interaction matrix for the grain boundary 7? is defined by 

- "J^AB + f^BA ■ 


The coefficients of the matrix R in Eq. (26) are obtained from R. as follows: 





0 

1 


Thus, for the example problem, E^ is given by 

= i|rX-rXf+ 0|r2N2-rMr 
i|rX-rXf+ i|rX-rXr- 


6. Extension of the variational problem to include grain boundaries 

Variational formulations of the Gurtin (2002) model for strain-gradient, single-crystal plas¬ 
ticity are investigated in Gurtin and Reddy (2014); Reddy (2011b) (see Reddy, 2011a, for the 
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Figure 3: Illustrative example of two grains 'Va and separated by a grain boundary 



polycrystalline case). The rate-dependent and independent cases are considered. The variational 
formulation developed in Reddy (2011b) is extended here to include the grain boundary. The 
associated incremental minimization problem, discussed in Sec. 6.3, provides the point of depar¬ 
ture for the numerical implementation in Sec. 7. 


6.1. The variational problem 

The spaces of displacement V and slips Q are defined by 

y = |h : Ui, e L^CV), 5m = 0 on , 

<2 = {r“ : ^ G L\y), = o on cTVh} • 

Note that the slips are continuous only in the grain and are, in general, discontinuous on the grain 
boundary. For the sake of simplicity, we assume micro-free conditions on d'Vf. 

The variational problem is obtained as follows. The strong form of the macroscopic equi¬ 
librium equation (14) and the microforce balance (15) are tested with an arbitrary displacement 
6u e V and slip 5y“ e Q and the result is integrated over the domain 'V. An integration by parts 
is performed leading to 


0 = I E(6u) : T dv - I Su ■ t* da, 


0 = J 6r‘^T"dv-J Sy^JT^dv- J V6y"-^dv + J 5r“[r ' «] da. 


(27) 

(28) 


The last term in Eq. (28), arising from the integration by parts, is the weak form of the microforce 

20 





balance on the grain boundary (17) evaluated on the interior boundaries of the grain, can be 
expressed as 


J ■«] dfl = J ■ «a ]+■ nsw da 

G G 

= f da 

G 

" ~f dfl 

G 




srK, 




‘B.disJ 


dfl. 


Using the above relation, the constitutive relations for the stress T (9) and the (energetic) 
vectorial microforce (10), the definition of the resolved shear stress, and the plastic flow relations 
in the bulk (22) and on the grain boundary (24) yields the variational problem: given the surface 
traction /*, find the displacement u e V and the plastic slips y“ € Q that satisfy 



for all Su eV and Sy“ e Q, where E^(u,y) - E(u) - E'‘’(y), p“ = p“(V 7 “) and G - G(yA,yB)- 


6.2. The incremental problem 

The time interval of interest 0 < f < T is partitioned into N subintervals as 0 = fo < 0 < 
■ ■ ■ < t]^ - T, with At - t„ - f„_i = T/N. The value of a quantity w at time t„ is denoted w„. 
The rate of change of a quantity is approximated using an Euler-backward difference scheme as 
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vv ss Aw/Af. The incremental problem is obtained by evaluating relations (29) and (30) at t„ as 




g\pe 

0 = 


dv 




^»J/e 

»= 


: S“ dv ■ 


- 6u ■ t* da, 


f r 

- civ- 

J >57“ ^ J 

T 'V 


dv- Vdr^-^^dv 


"V 

5Dvis 


^ _ _ ^vp r' _ _ 

(5r“ [N“ - N“] : ^ dfl - dr“ [N“ - N“] : 

^ dG n J 


dG 


dfl, 


(31) 


(32) 


where 


dD“ 


dy‘ 

5Dvis 


Ay^ 




dG 




Atro 

|AG| 


At Gq 


sgnAy“, 

AG 

|AG| ’ 


and AG can be expressed as 


AG = ;^[Ar“N“-Ar“N“] . 


6.3. The minimization problem 

The functional J{u, y) is defined by 


J{u 


,r):= JT'^(£‘=) dv + J'Pd(p)dv + Af^ JZ)“^(Ar“/Af) dv 

'V V “ ^ 

+ J T(G)dfl + Af J Dvis(AG/Af) dfl. 


r 

— t* ‘ u 

n 


dfl 


a^iv 


(33) 


Proposition. The solution to the incremental problem (31)-(32) is a minimizer of the functional 
J\ that is, if w„ (u„, y„) is kinematically admissible and satisfies 


J(Un,yn) < J(u„,y„) (34) 

for all w ;= (un,y„) eVxQ, then (m„, 7 „) is a solution of the incremental problem. 


Proof. For convenience the subscript n is dropped. An arbitrary but admissible G is denoted by 
G = G(yA,yB)- Since w is a minimizer of J it follows from Eq. (34) that 

J(\n) < J([l - 0]\N + 0w) (35) 
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for arbitrary but admissible w with 0 < 0 < 1. Set 

f(AyVAt) = J D"JAyVAt) dv and J(AG/At) = J D^AG/At) da . 

"V G 

From the definition (34) of the functional J and the convexity of and D^i^, Eq. (35) 
becomes 

I 'P (w) dv + Af ^ f (Ay‘^ I At) + At] [AG/At) 

. - <,w . ,W) d,.. - OIA,E / (Ar"/A,) . OA,2 (A^/AO 






-0 j t* ■\u-u]da + [l-0]Atj[AG/At) + 0Atj[AG/At) . 

tyvN 

Rearrangement of the terms in Eq. (36) leads to the inequality 

n ['I' ([1 - 0] w + 0 w) - 'P (w)] dv + Af ^ [/ (AY/At) - f (Ar7A0] 

- J t* ■ \u - u] da + At[j (AG/At) - J (aG/A?)] > 0. 


trvN 


Noting that 


trvN 


w - W = [W - W„_i] - [W - W„_i] 


Aw 


Aw = w, 


Eq. (38) becomes 

J" ^ : [w - Aw] dv + Af ^ [f (Y/^t) - f (Ay“/A0] - J 

ry O' 

+ Af[7(G/Af)-7(AG/Af)] >0. 


t* ■\u - Ah] dfl 


d^N 


(36) 


(37) 


By letting 0 go to 0 and using the definition of the derivative, Eq. (37) becomes 
r ^ ; [W - w] dv + Af ^ [f (A7“/A0 - f (Ay“/Af)] 

-J f* -[h-h] dfl + At (AG/At)(AG/Af)] >0. (38) 


(39) 
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Now 


5^ ^ S'?" 5'1'd 5/ _ 

5w dE^^ ^ dpP dVY' 

= r ; ^ ■ Vt^ . 

a 


Hence, using the additive decomposition of E, Eq. (39) becomes 
0< Jr: £(m) - £(Am) - ^ [7" - Ay“] S“ dv + J ^ ■ V(^ - Ay“) dv 






• Af ^ [/ (T/^O - f (Ar"/A0] - r <* ■ [S'- Ah] dfl + Af [; (C/Af) - j (AG/Af)] . 


(40) 


Now set ^ s Ay“, which implies G = AG, to get 


J T : [E(u) - ^(Ah)] dv - J t* ■ [u - Ah] da > 0. 


Choose H = +Ah to get 


J T : E(u) dv = J /* ■ H da, 

/TVn 


"V 


which is (31) in the incremental problem. This leaves in Eq. (40) the inequality 
- rr:S“;^[r-Ar“] dv+ r^r-v6^-Ar“) dv 

fy O' 

ry ^ ry ^ 

+ Af 2 [/ (Y/At) - f (AY/m + At \j (G/Af) - J (AG/Af)] > 0. 


Choose^ = [1 - 0]Ay" + GY with 0 < 6 < 1 to obtain 




0 [Y - Ar“] dv ■ 


/? 


r-V(0[r-Ar“]) dv 


'V 


'V 


+ Af ^ [/ ([1 - 0]Ar“/Af - GY/At) - f (Ay"/At)] 

a 

+ At [] ([1 - 0]AG/Af - GG/At) - ] (AG/Af)] > 0. 
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Dividing through by 6 and letting 0 go to zero 


- r^r^r-Ar"] dv+ r 

ry a ry ^ 


O' ^ At OCj At 


dJ\ 


Choose ^ + Ay“ 


The final two terms can be expanded as 


At 0(j 


G = 0. 


y ^ 

Z—i Qja 


= V r 


80“. 


dG 


Ay" I (5y‘ 

A? a ^ ^ 


Ay*’ 


dv, 


"V 

5D„ 


5G 


: Gdfl. 


(41) 


Finally, substituting this result into Eq. (41) gives 

J ^r'^ydv- f J ;r‘Tdv-J^ J ^ : [N“ - N“]y da = 0. 


□ 


Remark. The variational problem and corresponding minimization formulation are easily modi¬ 
fied to account for the alternative form for the dissipation and flow relation presented in Sections 
5.3 - 5.5. The only terms affected are those involving integrals over the grain boundary; thus, the 
last term in (30), repeated here for convenience, is 


- f 

G 



K 


da , 


and this term becomes, using the definitions in Section 5.4, 



- dfl. 


Likewise, in the functional (33) for the minimization problem, the term 




-i-Af Dvis(AG/Af) dfl. 
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is replaced by 


+Af 




£)(AG“^/AO 


da. 


7. Finite element approximation 

The spatially discrete form of the minimization problem (33) is solved approximately us¬ 
ing the finite element method in conjunction with a Newton-Raphson procedure. The software 
AceGen (Korelc, 2002) is used to describe the finite element interpolation, and to compute the 
residual and (algorithmically consistent) tangent contributions directly from the prescribed func¬ 
tional (33) using automatic differentiation, at the level of the quadrature point. This approach 
ensures quadratic convergence of the algorithm and greatly simplifies the implementation. 

The bulk material is discretized using 8-noded hexahedral elements. The nodal degrees of 
freedom are (m, y). Elements containing a grain boundary are discretized using the same element 
as in the bulk, but with the dimension in the direction of the grain-boundary normal far less than 
the other dimensions. Importantly however, all quantities on the grain boundary are computed 
and integration is performed with respect to the mid-plane of the grain-boundary element. Thus 
the grain boundary is treated as a two dimensional manifold embedded in three-dimensional 
space. The kinematic coherence constraint (4) is enforced by adding the following penalty term 
to the functional J\ 




where k'^° > 0 is the penalty parameter. An alternative approach would be to double the number 
of slip degrees of freedom at the face of an element that forms the grain boundary. Given that a 
grain boundary is only 1-2 atoms thick it would be preferable to use this approach. However, the 
approach is inconvenient to implement within the particular finite element framework used here. 

The finite thickness of the grain boundary has a negligible influence on the numerical re¬ 
sults. The contributions from an element containing a grain boundary to the balance of linear 
momentum in the bulk is ignored. This is a reasonable assumption provided the thickness of 
the grain boundary element is chosen to be sufficiently small. Additionally an extension to ac¬ 
count for possible opening and sliding (see e.g. Gurtin and Anand, 2008) is possible within this 
framework. 

8. Numerical examples 

Three example problems are presented to illustrate the key features of the grain-boundary 
model. The first example illustrates the influence of misorientation of the crystal lattice between 
adjacent grains in a bi-crystal. The second investigates the influence of the orientation of the 
grain boundary in a bi-crystal. The third example is that of a face-centered-cubic polycrystal 
subject to tensile loading. 

The constitutive parameters, chosen where possible to match those in van Beers et al. (2013), 
are summarised in Table 1. The thickness of the grain-boundary elements is fixed at le-5 [xm. 
Energetic contributions at the grain boundary are not accounted for in the numerical examples. 
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Table 1: Constitutive parameters used for the numerical examples unless stated otherwise. When a parameter has been 
varied, the alternative value is given in braces. 


First Lame parameter 

A 

1.05 X 10-' 

N/pm 2 

Second Lame parameter 


5.4 X 10-2 

N/pm 2 

Burgers vector length 

b 

2.5 X lO-'' 

fim 

Back stress cut-off radius 

R 

2x 10' 

fim 

Reference slip rate 

to 

1 

/S 

Go 

1 

/s 

Slip resistance 

5 

1 X 10-2 

N/pm 2 


1 X 10-' (1 X 10') 

S 

N/pm 2 

Rate sensitivity 

m, in 

1 


Penalty parameter 


1 X 102 



8.1. Misorientation of crystal lattice between adjacent grains 

Consider the [50]^ i^im bi-crystal subject to shear-type loading shown in Fig. 4. Single slip 
is assumed in both grains A (lower) and B (upper). The crystal lattice in grain A is fixed at 
{s, m) - (^ 1 ,^ 2 ). The normal to the grain boundary is fixed at « = 62 - The crystal lattice in grain 
B is initially chosen as {s, m) - { 61 , 62 ) and then rotated in increments through 90° about the 63 
axis. Note, the problem is similar to the one discussed in Sec. 2.3. The size of the bi-crystal 
is small relative to the internal length scale of the gradient-plasticity model governing the bulk. 
Gradient terms will thus play a role. 



Figure 4: The domain of the misorientation problem and the computational mesh. 

Each grain is discretized using 10 x 40 x 8 elements. The lower face of grain A is prevented 
from displacing in the 62 direction. The lines formed from the intersection of the lower face and 
the left and, separately, the back face are prevented from displacing in the 61 and 63 directions, 
respectively. The upper face is translated 5 pm in the Cj direction. In addition, micro-hard 
boundary conditions, required for the well-posedness of the variational problem, are imposed on 
the upper and lower faces. 
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The various interaction moduli for the range of orientations of slip system B, shown in 
Fig. 5 (e), can be extracted from Fig. 2 by setting ag - 0. The interaction moduli 
Cab are zero as tha x n = 0. The intra-action modulus Cbb varies smoothly between 0 and 1 as 
ub is varied. 

The variation of y over the vertical line jA-ys passing through the centre of the bi-crystal yg 
for Ob - {0°, 15°, 30°, 45°, 60°, 75°, 90°) corresponding the micro-hard, Gurtin I and II, and 
the micro-free grain-boundary models is shown in Fig. 5. Due to the relatively small specimen 
size, a significant boundary layer is observed in the upper grain for all ub + 0. For ub - 0, 
no dislocation pile-up occurs at the grain boundary or the external boundary as the slip system 
in grain B is parallel to these surfaces. The slip system in grain A is fixed parallel to the lower 
surface and the grain boundary. A small boundary layer is thus present in the vicinity of the 
lower surface. 

The interaction moduli Caa = Cab = 0 result in similar distributions of slip in grain A for the 
Gurtin I and micro-free grain-boundary models. The micro-hard and Gurtin II models produce 
a near identical distribution of slip in grain A. The two Gurtin models for the slip in grain A 
differ only in the vicinity of the grain boundary as the Gurtin II model does not account for the 
interaction of the slip system and the grain boundary directly in the flow relation. 

The value of the slip resistance on the grain boundary controls, in part, the flux of dislocations 
across the grain boundary. The default value of 5 = 1 x 10 ' is relatively low and hence the 
flux of dislocations into grain B is relatively high. Selecting a higher value of 5 = 1 x lO', 
significantly decreases the flux of dislocations into the grain boundary predicted by the two 
Gurtin models, as shown in Fig. 6. 

The micro-free grain boundary model is imposed by setting the grain boundary slip resistance 
S - 0. The micro-force balance on the grain boundary (17) thus reduces to the homogeneous 
Neumann condition. 

In summary, the Gurtin II model influences the flux of dislocations only through the grain 
boundary slip resistance. The Gurtin I model accounts for this and the interaction between the 
adjacent grains and the grain boundary. This observation suggests that the grain boundary slip 
resistance in the Gurtin II model could be chosen as a function of the mismatch in slip system 
orientation at the grain boundary, thereby allowing high-angle grain boundaries to act as micro- 
hard. 

8.2. Orientation of the grain boundary 

The influence of the orientation of the grain-boundary is investigated by fixing the crystal 
lattice in both grains of the bi-crystal considered in Sec. 8.1 and varying the grain-boundary 
orientation ag between -25° and 25°, see Fig. 7. In addition, the contribution of the gradient 
terms in the bulk material is investigated by considering two bi-crystal sizes: a large [1]^ = [100]^ 
pm and a smaller [I]^ = [50]^ pm. The imposed displacement on the upper surface is 1/10. The 
slip direction and slip plane normal for grains A and B are obtained by rotating the basis vectors 
ei and 62 by -5° and 10°, respectively. The variation in the inter- and intra-grain interaction 
moduli for the range of ag is given in Fig. 7. The angle between the normals to the slip systems 
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Figure 5: Distribution of y over the line y'A-y'B for various choices of ffg in the misorientation problem for the different 
grain-boundary models, (a)-(d), with S = 1 X 10“*. The interaction moduli for the range of orientations of slip system 
B is shown in (e). 


and the grain boundary is relatively small for the range of ag used, resulting in relatively low 
values for the interaction moduli. The bi-crystal is discretized using 10 x 40 x 8 elements. 

The influence of specimen size for the various grain-boundary models shown in Fig. 8. The 
boundary layers increase with decreasing size. The growth of the boundary layer is more marked 
for the micro-hard model. 

For ag - 10°, only ^ 0, and the micro-free and Gurtin I models produce near-identical 
results. As in the previous example, the relatively small value for S - 1 x 10 ' results in similar 
distributions of y for the micro-free and Gurtin I models. The Gurtin II model behaves more like 
the micro-hard model. 

8.3. Multiple slip in a polycrystal 

Consider a [100]^ p,m polycrystal subject to tensile loading. The polycrystal is composed 
of 27 equal sized grains. The crystal lattice in each grain is that of a face-centered-cubic ma¬ 
terial with the orientation of the 12 slip systems relative to the unit cell shown in Table 2. The 
orientation of the unit cell in each grain is listed in Table 3. 

The purpose of the polycrystal example is to explore the various choices of grain bound¬ 
ary models for reasonably complex, three-dimensional problems. As seen for the single-slip 
examples, a micro-free grain boundary model allows for a high flux of dislocations while the 
micro-hard model acts as an impenetrable barrier resulting in dislocation pile-up. The ability 
of the two Gurtin grain boundary models to impede dislocation flow is controlled, in part, by 
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Figure 6: Distribution of y over the line y’A-y'B for various choices of ffg in the misorientation problem for the different 
grain-boundary models, (a)-(d), with S = 1 X 10*. 


the slip resistance of the grain boundary. The Gurtin I model also accounts for the slip system 
interaction at the grain boundary. The ability of the two Gurtin models to capture the spectrum of 
behaviour between the micro-hard and micro-free limits is investigated by examining a range of 
grain boundary resistances from a low 5 = 1 x 10 to an extremely high value of 5 = 1 x 10^°. 
The alternative measure of the mismatch at the grain boundary developed in Sec. 5.3 is also 
examined. 

The discretization of the domain and the boundary conditions are shown in Fig. 9. Each 
grain is discretized with 6^ elements. A displacement of = 5 ^im is applied on the boundary 
with outward normal n - [0,1,0]. The opposite boundary with outward normal n - [0,-1,0] 
is prevented from displacing in the y-direction. The additional constraints to prevent rigid body 
motion are indicated in Fig. 9. 

The y-component of the resultant traction on the right boundary for the Gurtin I, micro-hard 
and micro-free models is plotted against the prescribed displacement in Fig. 10 (a). As expected 
the micro-free condition provides a lower bound for the Gurtin models for very low values of 
S. From the theory and the numerical investigations involving single slip it is reasonable to 
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Figure 7: The domain of the grain-boundary orientation problem and the computational mesh. Also shown is the variation 
in the inter- and intra-action moduli for a range of ag. 


micro-hard I Gurtin II Gurtin I I I micro-free 

|l00/im| 





Figure 8: Distribution of y over the line yA-ys for various choices of ag and specimen size in the grain-boundary 
orientation problem. 
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Table 2: The orientation of the slip planes relative to the unit cell. 


a 



a 


nf 

a 



1 

[110] 

(111) 

5 

[10 1] 

(III) 

9 

[Oil] 

(III) 

2 

[10 1] 

(111) 

6 

[Oil] 

(111) 

10 

[110] 

(111) 

3 

[Oil] 

(111) 

7 

[110] 

(111) 

11 

[10 1] 

(111) 

4 

[110] 

(111) 

8 

[10 1] 

(111) 

12 

[0 11] 

(111) 


Table 3: The Euler angles, expressed as <!>, ^ 2 ] following the convention of Bunge (1969), for each of the grains 

in the polycrystal. 


Grain 

<An 

Grain 


Grain 

-AH 

1 

[243,166,220] 

10 

[129,160,356] 

19 

[0,154,252] 

2 

[359,171,65] 

11 

[200,103,71] 

20 

[239,78,296] 

3 

[260,142,153] 

12 

[29,35,196] 

21 

[123,139,0] 

4 

[146,32,131] 

13 

[87,15,331] 

22 

[0,85,106] 

5 

[258,7,296] 

14 

[337,133,240] 

23 

[39,166,311] 

6 

[100,133,240] 

15 

[315,38,36] 

24 

[40,134,327] 

7 

[328,74,214] 

16 

[97,161,67] 

25 

[135,158,14] 

8 

[65,46,98] 

17 

[178,39,276] 

26 

[98,99,258] 

9 

[160,86,359] 

18 

[250,152,211] 

27 

[275,79,118] 


expect the micro-hard condition to be an upper bound that is approached with increasing S. It’s 
important to note that the extremely high values of S chosen are not physically motivated, rather 
they penalize the response at the grain boundary. For the high-angle grain boundaries present 
in the current example it is physically reasonable that the models be capable of producing a 
response close to micro-hard. It is clear that this is not the case for the Gurtin I model. 

Further insight into the behaviour of the Gurtin I model can be obtained from the distribution 
of the dissipation Dyis, and the rate of slip on each of the grain boundaries for a high value of 5 = 
1 X lO'®, as shown in Fig. 11. As discussed in Sec. 5.2, the Gurtin I model allows for a recombi¬ 
nation of slip via the definition of the grain boundary Burgers tensor used to parametrize the flow 
relation. The numerical solution procedure discussed in Sec. 6 is based on the minimization of 
an incremental potential. The Gurtin I model therefore permits a solution where Dyis is approxi¬ 
mately zero on the grain boundary for non-zero slip rates, as seen in Fig. 11. The recombination 
of slip at the grain boundary in the Gurtin I model is the reason that the micro-hard response can 
not be achieved. It should be emphasised that this is not necessarily a limitation, but rather a 
feature of the model. This feature in the original formulation is that for non-zero slip on the grain 
boundary, the dissipation can be zero. Although not shown here, the same behaviour is exhibited 
when energetic contributions on the grain boundary are included. 
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Figure 9: Discretization of the polycrsytal composed of 27 grains. The numbering corresponds to the grains listed in 
Table 3. 


The Guitin II model can capture the range of responses from micro-free to micro-hard, 
as shown in Fig. 10 (c). The same is true for the modified grain boundary model shown in 
Fig. 10 (b). 
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Figure 10: Applied displacement versus the >»-component of resultant force on the right boundary for various grain¬ 
boundary models. The results of the Guitin I model are shown in (a), the modified formulation in (b), and the Gurtin II 
model in (c). 
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9. Discussion and conclusion 


The grain-boundary theory of Gurtin (2008b) has been compared to recently developed mod¬ 
els in the literature. The theories are equivalent for planar problems. The weak form of the 
governing equations has been reformulated as an incremental minimization problem. This re¬ 
formulation provides an efficient and elegant basis for the numerical implementation. A series 
of three-dimensional numerical examples elucidated the theory for single slip in a bi-crystal. 
Various features of the Gurtin (2008b) grain-boundary models were illustrated for a polycrystal 
composed of grains with a face-centered-cubic structure. 

The Gurtin I model captures the geometric complexity of the grain boundary. A feature of 
the model is that it does not capture the full range of responses between micro-hard and micro- 
free. That is, it does not reproduce the widely used micro-hard limit when the grain-boundary 
slip resistance is used to penalize dislocation flow. The Gurtin II model does not contain the 
geometrical information concerning the grain boundary but can reproduce the micro-hard limit. 
The alternative formulation based on a modified definition of the grain boundary Burgers ten¬ 
sor is capable of capturing the micro-hard response and contains geometrical information. The 
modified formulation reduces to the Gurtin I model for single slip. 

The computational efficiency of the formulation is greatly impacted by the choice of the 
grain boundary flow relation. The scaled calculation times relative to the micro-hard condition 
for the polycrystal example are: 1.7, 4.4, 24 (for the Gurtin II, Gurtin I, and modified models 
respectively). The implementation of the modified model needs to be optimized, but the cost 
appears prohibitive. 

The numerical simulations provide valuable insight into the models. They do not, however, 
allow one to judge the physical correctness of the model. A key challenge is therefore the vali¬ 
dation and calibration of this and other grain-boundary models using well-devised experiments 
and microscopic modelling approaches (e.g. dislocation dynamics). 

The extension of the Gurtin (2008b) theory to the finite-strain regime is the subject of a 
companion paper in preparation. One significant challenge is the construction of a finite-strain 
counterpart to the grain-boundary Burgers tensor. The bulk theory is well understood (Gurtin, 
2008a). 

Elastic effects at the grain boundary have been added to the model presented here by merging 
it with the surface elasticity theory of Gurtin and Murdoch (1975). Results will be presented in a 
separate contribution. The inclusion of thermal effects into the grain-boundary gradient-plasticity 
model would be a challenging and recommended extension. One possible extension has recently 
been considered by Bargmann and Ekh (2013). 
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Appendix A. Gurtin I model for model slip 


The following example demonstrates how the Gurtin I model acts for multi-slip. Recall that 
the dissipation function and the rate of change of the grain-boundary Burgers tensor are given by 

-im+l 

3 ., 

GoJ 

a 

Assume the double-slip problem shown in Fig. A. 12. 




Figure A. 12: Double slip example. 


The problem is planar. That is, all slip plane normals m^,m^,m^ and m|, slip directions 
s^, s^, and s^ and the grain boundary normal n are in the x-y-plane. The lattice of grain B is 

the lattice of grain A rotated around the z-axis. One component of G is the projection tensor 
that is 

_j [ 1 

Na = ® (iWa X n) = ® 

All entries in the z-direction are zero, giving 


0 0 

0 0 0 

The term s\^(fn\jiy - m}y^Ti^) is replaced with a and s\y(m^y^Tiy - with b and the same 
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—a 

replacement is done for the other three N;. Then the rate of the Burger tensor can be written as 



O O 

O O 

+fl 

0 0 c 

0 0 d 

+ 75 

0 0 e 

0 0 / 

+ 71 

o o 

o o 


o 

o 

o 


0 0 0 


0 0 0 


o 

o 

o 


It is clear that there are combinations of non-zero 7 “ leading to zero dissipation at the grain 
boundary. 
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